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' Abstract. With increasing energy the diamagnetic hydrogen atom undergoes a 

transition from regular to chaotic classical dynamics, and the closed orbits pass 
through various cascades of bifurcations. Closed orbit theory allows for the 
semiclassical calculation of photoabsorption spectra of the diamagnetic hydrogen atom. 
However, at the bifurcations the closed orbit contributions diverge. The singularities 
can be removed with the help of uniform semiclassical approximations which are 
^ , constructed over a wide energy range for different types of codimension one and two 

I catastrophes. Using the uniform approximations and applying the high-resolution 

■ harmonic inversion method we calculate fully resolved semiclassical photoabsorption 

, spectra, i.e., individual eigenenergies and transition matrix elements at laboratory 

magnetic field strengths, and compare them with the results of exact quantum 
calculations. 



o 
o 
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^ ' 1. Introduction 



Rydberg atoms in a magnetic field have become prototype examples of a quantum 
system with an underlying classical dynamics changing from regular to chaotic motion 
with increasing excitation energy [1-3]. The Garton-Tomkins resonances originally 
found in barium atoms [4] and similar types of experimentally observed long-range 
modulations [5-8] can be associated with classical closed orbits starting at and returning 
back to the nucleus. A deeper quantitative analysis and interpretation of these features 
is possible within semiclassical theories, such as periodic orbit theory [9] and, as a variant 
for the photoabsorption of atomic systems, closed orbit theory [10,11]. In these theories, 
either the density of states or atomic photoabsorption spectra, as functions of energy, are 
given as superpositions of a smoothly varying part and sinusoidal modulations, whose 
frequencies, amplitudes, and phases are given in terms of the classical parameters of the 
orbits. 
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While closed orbit theory has been successfully applied to the interpretation of 
quantum spectra in terms of the closed orbits of the underlying classical system [12-15], 
the inverse procedure, i.e., the semiclassical calculation of the energies and transition 
strengths of individual eigenstates, is much more challenging for the following reasons: 
Firstly, closed orbit theory formally requires the knowledge of the infinite set of all closed 
orbits, which is impossible to obtain for nonintegrable systems where the orbits must be 
searched numerically. Secondly, in both periodic orbit and closed orbit theory the infinite 
sum over all orbit contributions suffers from fundamental convergence problems, and, 
thirdly, in generic systems the orbits undergo bifurcations when the energy is varied, and 
the contributions of isolated orbits exhibit unphysical singularities at these bifurcation 
points. 

Substantial progress has already been achieved to overcome these problems 
separately: On the one hand, the bifurcations of closed and periodic orbits of the 
diamagnetic hydrogen atom have been investigated in [16-18] and the divergences 
of isolated orbit contributions at bifurcations have been removed with the help of 
uniform semiclassical approximations for various types of bifurcations of codimension 
one and two [19-22]. On the other hand, the harmonic inversion technique based on 
high-resolution signal processing has been introduced as a method for semiclassical 
quantization of generic systems [23-25]. This method allows one to extract discrete 
eigcncnergies and matrix elements from a finite set of classical orbits, and thereby 
circumvents the convergence problems of the infinite sums in closed orbit or periodic 
orbit theory. In its original form it was applied to spectra of the diamagnetic hydrogen 
atom at constant scaled energy [26, 27] where the need to account for the effects 
of bifurcations does not arise. A semiclassical quantization with bifurcating orbits 
that merges these independent strands of research has only recently been achieved 
for an integrable atomic system, viz. the hydrogen atom in an electric field [28,29]. 
Energy dependent photoabsorption Stark spectra have been obtained by considering 
the bifurcations of the "uphill" and "downhill" orbit parallel and antiparallel to the 
direction of the external electric field. 

In this paper we demonstrate that the semiclassical quantization with bifurcating 
orbits can be successfully applied to a more challenging system, viz. the hydrogen 
atom in a magnetic field, where the classical equations of motion are nonseparable and 
the dynamics undergoes a transition from regular to chaotic dynamics with increasing 
energy. The bifurcation scenarios encountered there are much more complicated than 
those in the hydrogen atom in an electric field, and different types of catastrophes with 
codimension one and two must be used to remove the divergences at the bifurcations. 

Although the numerical effort for the calculation of the semiclassical photoabsorp- 
tion spectrum of the diamagnetic hydrogen atom is much higher than for the corre- 
sponding exact quantum computations, the results of this paper are of fundamental 
interest for the development, understanding, and practical applications of semiclassi- 
cal theories. As Einstein [30] pointed out as early as 1917, the "old" quantum theory 
based on the Bohr-Sommerfeld quantization rules is doomed to failure when applied to 
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nonintegrable systems. About ninety years later we have now succeeded in obtaining 
the high-resolution photoabsorption spectra of a nontrivial atomic system with mixed 
regular-chaotic dynamics semiclassically from first principles. The necessary ingredients 
are closed orbit theory, uniform semiclassical approximations at bifurcations, and the 
harmonic inversion method. 

The paper is organized as follows. In Sec. El the classical dynamics of the hydrogen 
atom in a magnetic field and various types of closed orbit bifurcations are discussed. In 
Sec. 01 closed orbit theory is introduced and the uniform approximations at bifurcations 
of closed orbits are constructed. Semiclassical high-resolution photoabsorption spectra 
with individual eigenenergies and transition matrix elements are obtained by application 
of the harmonic inversion method in Sec. |3] and are compared with exact quantum 
spectra. Concluding remarks are given in Sec. IHl 



2. Classical dynamics and closed orbit bifurcations 

The classical dynamics of the diamagnetic Kepler problem has already been discussed 
extensively in the literature (for reviews see, e.g., [1-3]). Here we briefly recapitulate 
the essentials which are necessary to what follows. The Hamiltonian in atomic units 
(with the magnetic field B = Be^ along the 2;-axis, 7 = i?/(2.35 x 10^ T), and angular 
momentum component = 0) reads 

H = \v'-\ + \^\x' + y')=E, (1) 

where E is the energy. Using semiparabolic coordinates /x = \/r + z, v = ^r — z the 
Hamiltonian (0) can be transformed to 

h = lipl+pl) - Ei^^^ + u') + l^'f^'Af^' + u^)=2. (2) 

Note that Hamilton's equations of motion derived from Q are free of singularities at 
the Coulomb centre. The Hamiltonian (Q) is invariant under a reflection in the xy-plane 
perpendicular to the magnetic field. As a consequence, all orbits that are not in that 
plane occur in symmetric pairs. The closed orbits leave the nucleus (r = 0) with an 
initial angle to the ^-axis and return back to the origin with final angle after 
time period T. The stability properties of the closed orbits are given in terms of the 
2x2 monodromy matrix M, which linearly maps local deviations {6q, Sp) of the starting 
point in the directions perpendicular to the orbit in coordinate and momentum space 
onto local deviations of the final point: 

\6piT)J \5pi0)) \m2i m22j\6pWJ' 

Closed orbits bifurcate when the element of the monodromy matrix vanishes. (Note 
that a different condition det(M — 1) = is valid for periodic orbit bifurcations.) 

The classical dynamics does not depend on the energy E and magnetic field strength 
7 separately. Instead, the scaled Hamiltonian H = H'j~'^^^ is independent of 7 if it is 
expressed in terms of the scaled semiparabolic coordinates fl = 7^/^/i and u = j^^^u, 
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so that the scaled energy E = £'7"^/^ is the only control parameter. The time and 
classical action scale as t = ^7 and s = 57^^^, respectively, and the matrix element mi2 
of the monodromy matrix M in equation (jHl), which is important in closed orbit theory, 
scales as mi2 = mi27^^^. 

In the limit of infinitely negative energy [E — > —00) only two closed orbits parallel 
and perpendicular to the magnetic field axis exist. Their multiple repetitions are called 
basic vibrators and rotators R^, respectively [6,7], where the index /i gives the number 
of repetitions. When the energy is increased, the basic vibrators and rotators undergo 
cascades of bifurcations where new closed orbits and R'j^ are created in a systematic 
way (see Sees. I2.l1 and l2.2l below). These orbits can run through further bifurcations as 
discussed in Sec. 12.31 Furthermore, new closed orbits which are not directly related to 
the bifurcation tree of the basic vibrators and rotators are created "out of nowhere" by 
tangent bifurcations. With increasing energy a transition from nearly regular to chaotic 
phase space takes place, along with a rapid proliferation of closed orbits, and thus the 
semiclassical quantization with bifurcating orbits becomes more and more challenging 
with growing energy. 

2.1. Bifurcations of the parallel orbit 

With increasing energy each repetition V"^ of the parallel orbit undergoes an infinite 
sequence of bifurcations with an accumulation point at the field free ionization threshold 
E = 0. Individual bifurcations are counted by integer numbers u. The bifurcations of 
the parallel orbit are illustrated in figure where the periods of the basic vibrators 
are plotted as functions of the scaled energy E. At the bifurcation points new orbits 
called [6,7] are created. The bifurcation points are marked by crosses in figure [T^, 
and some of them are labelled by the symbols of the bifurcating closed orbits. The 
periods of the bifurcating orbits as functions of the scaled energy are not shown to 
keep the figure concise. From figure ^ it becomes evident that bifurcations occur rather 
frequently in energy, and the correct handling of bifurcations is of crucial importance 
for the semiclassical quantization of this system. 

The bifurcations of the parallel orbit of the diamagnetic hydrogen atom resemble 
those of the "uphill" and "downhill" orbit of the hydrogen atom in an electric field 
[28,29,31]. The system is rotationally symmetric around the parallel (axial) orbit, 
and a bundle of three-dimensional non-axial orbits collide with the axial orbit in a 
pitchfork bifurcation. The non-axial orbits are real orbits or complex "ghost" orbits 
in the complex continuation of the phase space when the value of E is above or below 
the bifurcation energy, respectively. As an example, the bifurcation scenario of the 
orbits V4 and is illustrated in figure |21 The difference AS/27T = {S^.^ — S^o-a)/2n 
between the scaled actions of the axial (ax) and the non-axial (non) orbit and the 
element of the scaled monodromy matrix are shown as functions of the scaled 
energy. The solid and dashed lines refer to the real and ghost orbits, respectively. The 
approximately linear dependence of fhi2 on the energy and the quadratic behaviour 
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Figure 1. Periods of (a) the basic vibrators Vf, and (b) the basic rotators i?^. The 
bifurcations are marked by crosses, some of them are labeUed by the symbols or i?^ 
of the bifurcating closed orbits. The periods of the bifurcating orbits are not shown to 
keep the figure concise. 
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Figure 2. Closed orbit parameters for the pitchfork bifurcation of the orbits V4 and 
Vl . Sohd and dashed hnes mark the real and ghost orbits, respectively. 



of AS" around the bifurcation point is typical of pitchfork bifurcations. For the ghost 
orbits the initial angle di (not shown in figure ^ is purely imaginary. The classical 
action and monodromy matrix, however, are real valued. Note that the closed orbit 
parameters of both the real and ghost orbits are required for the construction of the 
uniform semiclassical approximations in Sec. El 
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2.2. Bifurcations of the perpendicular orbit 

The multiple repetitions of the orbit perpendicular to the magnetic field, i.e., the basic 
rotators i?^ also undergo a sequence of bifurcations, where new closed orbits R"^ are 
created [6, 7]. Contrary to the bifurcations of the basic vibrators, the number of 
bifurcations of a rotator is finite. The integer v that identifies the bifurcations is limited 
by z/ < /i. An overview of the bifurcations of the perpendicular orbit is given in figure^, 
where the bifurcation points are marked by crosses, with some of them labelled by the 
symbols W^. 

The bifurcations of the basic rotators have been investigated in [16-19] and turn 
out to be much more subtle than those of the basic vibrators. The scenario is always a 
sequence of a tangent and a pitchfork bifurcation which occur at two nearby bifurcation 
energies. As an example we discuss the creation of the closed orbit R\ in a bifurcation 
out of the third repetition of the basic rotator i?3. The closed orbit parameters are 
presented in figure El which shows the energy dependence of the action difference AS/ 271 
with the scaled action of the orbit R3 taken as the reference action, and the element mi2 
of the scaled monodromy matrix. Real and complex "ghost" orbits are marked by solid 
and dashed lines, respectively. The real closed orbits R^ and Rl collide in a pitchfork 
bifurcation at scaled energy E = —0.48284. Below that bifurcation energy a ghost 
orbit in the complex phase space, albeit with real action and monodromy matrix, does 
exist. This orbit participates in a tangent bifurcation at scaled energy E = —0.48477. 
Unlike a conventional tangent bifurcation, where two real orbits are created out of 
ghost orbit predecessors [32] , the tangent bifurcation shown in figure El possesses the 
peculiar property that all participating orbits are complex ghosts. One ghost orbit 
with complex action and monodromy matrix and its complex conjugate companion 
bifurcate a.t E = —0.48477 into two genuinely different ghost orbits with real actions 
and monodromy matrices. The significance of ghost orbit bifurcations in semiclassical 
spectra has already been demonstrated in [21,22]. The closed orbit parameters, including 
those of the complex ghosts, are required for the construction of uniform semiclassical 
approximations in Sec. El 

A systematic survey of the bifurcations of the basic rotators reveals that the 
bifurcation scenario with the ghost orbit tangent bifurcation as shown in figure El is 
restricted to bifurcations of the basic rotators at scaled energies E < —0.418 [33]. 
At higher energies conventional tangent bifurcations are found instead of ghost orbit 
bifurcations. Examples are the creation of orbits R^ at E ^ —0.317 [19] or at 
E ^ -0.209. 

2.3. Bifurcations of orbits V/^ and R^^ 

After being created, the vibrators and rotators R^^ can themselves undergo further 
bifurcations with increasing energy. These are pitchfork bifurcations, where a pair of 
asymmetric orbits (with different initial and final angles, di 7^ df) separate from a 
central symmetric obit with di = df 01 di = it — i^f, which is real below and above 




Figure 3. Closed orbit parameters for the bifurcation scenario where the orbit Rl 
is created from R^. A sequence of a tangent bifurcation where only ghost orbits 
participate and a nearby pitchfork bifurcation is the typical scenario for all bifurcations 
of the perpendicular orbit at scaled energies E < —0.418. Real and ghost orbits are 
marked by solid and dashed lines, respectively. 



the bifurcation. On the two asymmetric orbits the electron follows the same trajectory 
in different directions, i.e., the final angle of one orbit is the initial angle of the other 
orbit and vice versa. Below the bifurcation point the asymmetric orbits become a pair 
of complex conjugate ghost orbits. Again the electron follows the same trajectory in 
different directions and therefore must have real actions and monodromy matrices. For 
illustration the closed orbit parameters of the first pitchfork bifurcation of the orbit Vj^ 
a.t E = —0.514 are presented in figure^ A typical feature of the pitchfork bifurcations 
is the nearly linear dependence of the monodromy matrix element rhu and the nearly 
quadratic behaviour of the action difference AS around the bifurcation energy. For 
pitchfork bifurcations of vibrators the action of the asymmetric orbits exceeds the 
action of the symmetric orbit. Pitchfork bifurcations of the rotators Wj^ show the 
opposite behaviour, i.e., the action of the symmetric orbit i?^ exceeds the action of 
the newly created asymmetric orbits. 

The bifurcations discussed so far are sufficient to describe the complete bifurcation 
tree of all closed orbits at scaled energies E < —0.5 and with recurrence times 
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Figure 4. Closed orbit parameters for the pitchfork bifurcation of the vibrator V^j^ at 
E — —0.514. The real asymmetric orbit is labelled V^l . The predecessor ghost orbit 
is marked by dashed lines. 

T jln < 12. However, at higher energies more types of bifurcations exist, e.g., new 
closed orbits can be created in tangent bifurcations (without an accompanying pitchfork 
bifurcation as for the bifurcations of the perpendicular orbit discussed in Sec. I2.2|l . 
The simplest of these orbits, called Xi in [6,7,19], is created at the scaled energy 
E = —0.11544216. As the calculations of the semiclassical photoabsorption spectra in 
Sec. Hlare restricted to energies E < —0.5 we need not discuss the tangent bifurcations 
and the corresponding uniform semiclassical approximations of such orbits in more 
detail. 
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3. Closed orbit theory and uniform approximations 

Closed orbit theory [10, 11] provides a semiclassical approximation to photoabsorption 
spectra of atoms in external fields, where the electron is excited from a low lying initial 
state |\E'i) to a final Rydberg state \^n)- It is convenient to introduce the quantum 
mechanical response function 

g^-^(E) = --(^.IDG+DI^,) = --y (4) 

n 

with D the dipole operator. En the eigenenergy of the eigenstate \^n)^ and the 
retarded Green's function. From the response function (j3)) relevant physical data such 
as the oscillator strength 

f{E) = 2{E-E{)\mg^^{E) (5) 

are readily obtained. The semiclassical approximation to the exact quantum response 
function (jH) is given by closed orbit theory as a sum of a smooth background term and 
an oscillatory part 

= ^ Ao(^)e'(^^°(^)-5^^°) , (6) 

CO 

where the sum is to be taken over all closed orbits (co) starting at and returning back 
to the nucleus. S^o and /ico are the classical action and Maslov index of the closed orbit, 
respectively. The amplitudes depend on the symmetry of the orbits and read 



for non-axial closed orbits, which in the three-dimensional coordinate space form a 
rotationally invariant family of orbits around the field axis, and 

^^o = j^ymyw (s) 

|"^12| 

for axial orbits [12,34], i.e., the basic vibrators parallel to the magnetic field. In Eqs. (|7j) 
and (jHI) and i)f are the initial and final angle, and is an element of the monodromy 
matrix of the closed orbit. Note that these parameters depend on the energy E. The 
functions y{'&) characterize the initial state |\E'j) and the polarization of the dipole 
transition, and are linear combinations of spherical harmonics Yim{'&,0). The closed 
orbit amplitudes (|7j) and (jH)) are valid in integrable as well as chaotic regimes. By 
contrast, the trace formulas of periodic orbit theory are different for integrable [35] 
and chaotic [9] systems. However, the closed orbits must be isolated, i.e., Eqs. ((Tj) and 
(jHI) fail near bifurcations where different orbits approach each other and eventually 
collide. The element of the monodromy matrix vanishes at bifurcations, and 
the semiclassical amplitudes ^co of the isolated closed orbit contributions suffer from 
unphysical singularities. 

To obtain a smooth contribution to the semiclassical response function g°^'^{E) 
where closed orbits bifurcate, uniform approximations are needed. Their construction 
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requires a detailed description of the bifurcation scenario. In the language of catastrophe 
theory [36, 37], this can be achieved in terms of normal forms whose stationary 
points correspond to the classical closed orbits [19, 38]. The codimension of the 
bifurcation scenario coincides with the codimension of its normal form. For a generic 
bifurcation, it is at most the number of external parameters, which is one in the case 
of the diamagnetic Kepler problem. In this system two types of generic bifurcations 
exist: tangent bifurcations and pitchfork bifurcations. They can be described by the 
fold and the symmetric cusp, respectively. More complicated bifurcation scenarios 
are composed of several generic bifurcations and modelled by catastrophes of higher 
codimension. If the individual bifurcations are closely spaced, it is important to 
construct a uniform approximation that describes them collectively. Various uniform 
semiclassical approximations have already been constructed for the hydrogen atom in 
an electric [28,29,39,40] and a magnetic [19,20] field. 

As a starting point for catastrophes of corank one we consider the ansatz 

Ci(i^) = y"p(t)e'*-(*Mte^(^°-5^«), (9) 

where $a(^) is the normal form of the catastrophe depending on the parameters 
a = (ai,a2, . . . , ctfc) with k being the codimension. The uniform approximation is 
supposed to reproduce the closed orbit sum of all orbits participating in the bifurcation 
scenario if the distance from the bifurcation is large. 5*0 is an energy dependent reference 
action, e.g., the action of a central closed orbit. The integer z/q and the function p{t) 
must be chosen to asymptotically provide the correct phase and amplitude of the uniform 
approximation. The classical actions of the closed orbits contributing to the bifurcation 
scenario are related to the stationary values of the normal form $a(^)- To determine 
the parameters of the normal form $a(i) and the amplitude function p{t) we use the 
asymptotic expressions of the uniform approximation Q far away from the bifurcations, 
where the integral can be evaluated in stationary phase approximation. The stationary 
phase (sp) method applied to equation yields 

where the tn are the stationary points of $a(^)- The constant integer z/q is given by the 
Maslov indices of the closed orbits, which change at the bifurcations, i.e., /x^j = z/q + 
with 

I 0, otherwise. 

The sum over the stationary points in equation (fTUj) is identified with the sum © over 
the closed orbits colliding in the bifurcation scenario described by the normal form $a(^)- 
From the comparison of the two equations © and ITUI) we obtain the conditions 

So + ^(tn) = Sn (12) 
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and 

V27ri j9(t„) 



, An . (13) 

Eqs. (I12j) and ()13|1 are valid not only for real orbits but hold also for complex ghost 
orbits with the slight modifications [41] that \mi2\ and |$"(t„)| are replaced with 
sign(Remi2) mi2 and sign(Re $"(t„)) respectively. These equations are now used 

to determine, from the numerically calculated closed orbit parameters S'„ and An, the 
parameters of the normal form $a(^) and the amplitude functions p{t) for various types 
of catastrophes. The functions $a(^) and p{t) are then inserted in equation © to obtain 
the uniform approximations. In the following we discuss the butterfly catastrophe 
related to bifurcations of the perpendicular orbit, the symmetric cusp catastrophe 
related to bifurcations of the rotators and vibrators, and the uniform approximations 
for bifurcations of the parallel orbit. 



3.1. Uniform approximations for bifurcations of the perpendicular orbit 

The bifurcation scenarios of the perpendicular orbit discussed in Sec. 12.21 are described 
by the codimension-2 symmetric butterfly catastrophe [19]. A local approximation, 
which removes the singularities at the bifurcation energies but is not valid at energies 
far away from the bifurcation points has been presented in [19]. Here a uniform 
approximation is derived that removes the unphysical divergences at the bifurcations 
and agrees asymptotically with the closed orbit sum of the isolated orbits. 
The normal form of the symmetric butterfly catastrophe reads 

^^t) = -{t^ + xt^ + yt') , (14) 

where the two real unfolding parameters x and y depend on the energy and magnetic 
field strength and must be determined from equation (fT^ . To that end, the stationary 
points of $(t) at 

to = 0, ti,2 = ±V(-a; + 5)/3, t3,4 = - 5)/3 , (15) 

with 5 = a/x^ — 3?/ are identified with the (real or complex) closed orbits contributing 
to the bifurcation scenario: Real stationary points are real orbits, purely imaginary 
stationary points are ghost orbits with real action, and complex points are ghost orbits 
with complex action. The trivial stationary point to is the perpendicular orbit that is real 
on both sides of the bifurcation. The nontrivial stationary points are real, imaginary, or 
general complex numbers in various regions of the parameters x and y as illustrated in 
figureEl The lines y = and y = x^/3 divide the (x, ?/)-plane into different domains and 
characterize the parameter values where pitchfork or tangent bifurcations, respectively, 
of the closed orbits occur. The stationary points ti...4 are imaginary in region A, complex 
in region B and real in region C. In region D, ti 2 are real whereas ^3^4 are imaginary. 

These observations translate into the following bifurcation scenarios: In region D 
the stationary points ti^2 represent the real orbits R'^. If x > 0, as ?/ is increased 
these orbits collide with the perpendicular orbit (to) at y = and become imaginary. 
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Figure 5. Characterization of the stationary points of the normal form in 
equation l|14|) . The functions y = and y = divide the (x,y)-plane into different 
domains. A: ti..A G iK; B: ti..A € C; C: ti..A G R; D: ti^2 G K, ^3,4 G iM. 



corresponding to ghost orbits with real action, in region A. At |/ = they colhde with 
^3 4, and finally in region B the stationary points ^1. 4 represent a quadruple of ghost 
orbits with complex action. This scenario with unfolding parameter x > is observed 
for bifurcations at scaled energies E < —0.418. The scenario for a; < is similar, the 
difference being that all stationary points and corresponding closed orbits are real in 
region C. This scenario has been observed at scaled energies E > —0.418 [33]. The 
stationary values of the normal form in equation ()14j) and the second derivative 
$"(t) as functions of the unfolding parameter y are presented in figure IHl Evidently, 
figure El qualitatively agrees with the closed orbit parameters AS* and shown in 
figure El 

The unfolding parameters x and y can now be determined from the actions of 
the colhding orbits using equation (fT^ . When Sq is identified with the action of the 
perpendicular orbit, we obtain 

^ $1 = - A(x - 5)5' + ^xy = 5i - 5o = AS, , (16) 

$(t3) ^ $3 = + 5)^2 ^ = s,-So = AS, , (17) 

where 5*1 and S3 denote the actions of the corresponding orbits described above. The 
sum and the difference of Eqs. (fTHjl and (fTTjl yields 

ASi - AS3 = $1 - $3 = ^ (18) 
A^i + A53 = $1 + $3 = ^ (a:^ - 35'x) . (19) 
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Figure 6. (a) Stationary values of the normal form in equation 114|l as function of 
the unfolding parameter y (x = 3 > 0). (b) Second derivative at the stationary 

points. The qualitative agreement with figure 13 is evident. 



Now 



5= v/x2-3i/ = 3^(<l>i-<l>3)/4 (20) 

follows from equation ()18p. and the unfolding parameter x is obtained from the solution 
of the cubic equation (fTUj) using Cardano's formula as [42] 



($1 + $3) + 2V<l>i$3 



-4 



($1 + $3) - 2v/$i<f 3 ,(21) 



with A G {1, (—1 ± i-\/3)/2}. The parameter x becomes a smooth function of the energy 
with the choice 



A 



1, E<E,, 
;i + iy3)/2, E>E,, 



(22) 



where Ec is the bifurcation energy of the pitchfork bifurcation, i.e. $1 = 0. Finally the 
unfolding parameter y is given as ?/ = (x^ — 5^)/3. 

The next step is to construct the function p{t) in the uniform approximation 
in such a way that equation (fT^ is valid at the stationary points t„ of the normal 
form $(t). There is considerable freedom in doing so, and one will strive for an ansatz 
for p{t) that is as simple as possible. In all cases discussed in the literature so far 
(e.g., [38,40,41,43]), it has been found sufficiently accurate to choose pit) to be a 
low-order polynomial with as many undetermined coefficients as there are conditions 
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imposed by (fTS|) . For the bifurcation scenario described by (fT^ . however, we find that a 
polynomial ansatz for p{t) yields a uniform approximation that in the bifurcation region 
differs from the expected results and assumes the correct asymptotic behaviour only 
at huge distances from the bifurcation. This scenario thus calls for a more thorough 
analysis of the amplitude function. 

Because contributions to the integral in the uniform approximation Q arise chiefly 
in a neighbourhood of the origin, where the stationary points of the normal form are 
located, the amplitude function p{t) need only be accurate in that region. A polynomial 
ansatz is justifled if p{t) is nearly constant in the region of interest, so that the t- 
dependent terms are small correction to the constant term and high-order terms that are 
not included in the ansatz are negligible. A comparison of figure El and figure El reveals 
that the energy- dependence of the monodromy matrix element mi2 is well described 
by the second derivative of the normal form. To satisfy it remains for p{t) to 

describe the energy dependence of the amplitude through the angles 'di and 'dj. These 
considerations suggest the ansatz 

p{t) = {at'' + bt^ + c)y\m) . (23) 

with the same angular function y{'d) as in equation ((7|). The mapping {}{t) from the 
normal form coordinate t to the angle d is in turn modelled by the polynomial ansatz 



IT 



+ vt + uf 



(24) 



with the coefficients u and v chosen such that d{t) maps the stationary points of the 
normal form (fT^ onto the numerically determined angles di and df. (Note that di = df 
oi di = TT — d f for all orbits involved.) 

Inserting equation into (fT3|) . we obtain the three coefficients a, 6, c in as 




3x 
2^ 




37r 



Al^ 7]i5{x -5)- Al^7]z5{x + 5) 



A\^r]i6{x -6)+ Al^/r]i6{x + 5) 



:i-i) 



(25) 




Al^r]i5{x -5)+ Al^r]z5{x + 5) 
AUil^5{x-5)-Ayr]^6{x + 5)\ \ (1-i) , (26) 

(27) 

where ^* = An/ {y{di)y{d j)) denotes the semiclassical amplitudes in equation ((Zj) 
without the angular functions, x, y and 5 are the parameters as introduced above, 
and rjn = sign (Re $"(t,i)). Note that all coefficients in the normal form $(i:) and the 
amplitude function p{t) are now given as explicit functions of the closed orbit parameters 
of the real and complex (ghost) orbits involved in the bifurcation scenario. 




Figure 7. Absolute value of the semiclassical response function for the bifurcation 

scenario of orbits Rq and Dashed line: isolated closed orbit sum; solid line: uniform 
approximation. 



With the normal form and the amphtude function p{t) at hand, it is possible 
to evaluate the uniform approximation The integral must be solved numerically. 
For if: ±00 the integrand is highly oscillating and must therefore be regularized 
by multiplication with a factor of the form exp(— et™") with the small e > and the 
power m > chosen appropriately. As an example figure U\ presents the absolute 
value of the semiclassical response function for the orbits Rq, i?g, and the ghost orbits 
associated in the bifurcation scenario. The isolated closed orbit sum (dashed line) 
suffers from the unphysical divergence around E = —0.78. By contrast, the uniform 
approximation (solid line) is a smooth function at all energies. The modulations of 
the amplitude a.t E > —0.7 are caused by the interference of the real closed orbits -Re 
and i?g. Note that the uniform approximation at large distances from the bifurcation 
energies asymptotically agrees with the isolated closed orbit sum. 

3.2. Uniform approximations for bifurcations of the rotators and vibrators 

The pitchfork bifurcations of the rotators Rj^ and vibrators discussed in Sec. 12.31 are 
described by the normal form of the symmetrized cusp catastrophe 

$,(t) = \t' - \at^ , (28) 

which has one unfolding parameter a. The stationary points at t = and t = ±^/a 
correspond to the symmetric (sym) and asymmetric (asym) orbits, respectively. The 
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normal form parameter 
a = ±2^/5; 



sym 



asym 



(29) 



is a function of the actions of the closed orbits involved. It has to be chosen positive if 
the asymmetric orbits are real, and negative otherwise. 

Due to symmetry properties of the orbits the amplitude function p{t) in equation ^ 
must be an even function of t. It is sufficient to use a simple polynomial ansatz 
p(t) = Po + P2t'^. The coefficients are obtained as 



Po 



P2 




asym 



V2As 



sym 



(30) 
(31) 



with ^sym and ^asym the amplitudes of the (isolated) symmetric and asymmetric closed 
orbits as defined in equation ((7j). (Equations (jHUj) and (jHT|) slightly differ from formulae 
given in [38] due to a different handling of the Maslov phase in equations (jH)) and ((Tj).) 
The uniform approximation can now be written as 



osc 
i'uni 



p{t)e"^-^'Mte' 



(po/o + P2/2)e'(^-'"-5"°) 



(32) 



with the integrals Iq = j exp(z$a(t))dt and I2 = J exp(i$a(t))dt. The integrals can 
be evaluated analytically in terms of Bessel functions Ju{z) [42] and read [38] 



71 

2 



a\e 



\a\e 



1 

2a ~ 



1/4 



+ signae-*^/^Ji 



/4 



(33) 



1/4 



Sign a e 



+ signae-^''/Vi 



/4 



J-3/4 I — 



J5/4 I — 



(.34) 



The normal form (j28p and thus the uniform approximation (|32j) describe the pitchfork 
bifurcation when the action of the symmetric orbit is larger than the action of the 
asymmetric orbits in the vicinity of the bifurcation. This is true for the rotators i?^. 
For the vibrators the converse is true, i.e., close to the bifurcation the action of the 
asymmetric orbits exceeds the action of the symmetric orbit V^. In this case, called a 
dual cusp [36], the normal form $a(^) must be replaced with —^ait), which changes the 
sign of the stationary values. The uniform approximation for the dual cusp is obtained 
by replacing the integrals Iq and I2 in equation with its complex conjugate, i.e., 
Ci = {poI*o+P2l*2)exp{i{S,^^ - fz/o)). 



3.3. Uniform approximations for bifurcations of the parallel orbit 

In bifurcations of orbits parallel to the external field a rotationally symmetric bundle of 
non- axial orbits splits from the axial orbit. This is true for both an external magnetic 
and an electric field. The bifurcations of the axial orbits are described by the normal 
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form ^aif) = \t'^ — \0't'^ which is formally the symmetric cusp, equation (PH|). but 
with the difference that t is interpreted as a radial coordinate, t = a/x^ + y"^. For the 
amplitude function in the ansatz p(t) = po + Viif' ~ is used [29]. The stationary 
points at t = ^/a describing the bundle of non-axial orbits lie on a circle with centre at 
t = 0, which is the stationary point describing the isolated axial orbit. The parameter 
a = ±2\/S'ax — >S'non IS related to the classical action of the orbits, with positive and 
negative values of a referring to real and ghost orbits, respectively. 

The bifurcations of the "uphill" and "downhill" orbits parallel and antiparallel to 
the electric field axis in the Stark system have already been investigated and a uniform 
approximation for the creation or destruction of the non-axial orbits from the axial orbits 
has been constructed [39,40]. The uniform approximation is valid also to describe the 
bifurcations of the basic vibrators in a magnetic field, and can be written in the concise 
form [28, 29] 

(35) 



with 

1(a) = e-'"'/^ 

given in terms of the Fresnel integrals C{z) and S{z) [42]. (Similar as above equation (j35j) 
slightly differs from the result given in [28,29] due to a different handling of the Maslov 
phase in equations ©^(El)-) 



A 

r+i' 



-Ka) 



27r 




4. High resolution photoabsorption spectra 

With the closed orbit theory and the uniform approximations at hand we can now 
obtain the semiclassical response function g^^{E) via summation of the closed orbit 
contributions. In the vicinity of bifurcations the contributions of isolated orbits must be 
replaced with the uniform approximations. The semiclassical photoabsorption spectrum 
is then readily given by equation (0) with g'^^{E) replaced with its semiclassical 
analogue. However, the closed orbit sum diverges when it is extended over all closed 
orbits. When it is truncated, e.g., by neglecting orbits with recurrence time T > Tmax, 
it yields only low resolution spectra. To obtain high resolution spectra, i.e., discrete 
eigenenergies and individual transition matrix elements dn = \{^ i\D\^ ■, we 
adopt the harmonic inversion method [25] which has been successfully applied in 
semiclassical mechanics, either to extract the actions and amplitudes of classical orbits 
from quantum spectra [44] or to calculate quantum mechanical quantities from classical 
orbits [23,24,26,27]. It has also been demonstrated that harmonic inversion is a 
powerful tool for semiclassical quantization using bifurcating orbits [28,29]. To keep 
our presentation self-contained, we briefly outline the basic ideas. 
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The harmonic inversion method 

In a first step, both the quantum (jlj) and the semiclassical response function ^ - 
the smooth part can be neglected - are Fourier transformed into time domain. The 
Fourier integrals are restricted to the energy window [-Emin, -f'max] where the closed orbit 
parameters have been calculated. The windowed Fourier transforms result in the band- 
limited time signals 

^-(t) = / yAo(^)e'^^°(^)e-'^*dE . (38) 

ZTT Ip . 

■J ^mm CO 

In the quantum signal (jHTj) the sum is restricted to the eigenenergies in the range 
-Emin < En < -Emax, i-G., ouly a relatively small number of parameters {En,dn} must 
be determined if the energy window is chosen appropriately. In the semiclassical signal 
(I38|l only those closed orbits contribute within a stationary phase approximation whose 
recurrence times T are less than the total length T^ax of the time signal. This means 
that the semiclassical signal (jHHj) can be constructed if the set of the closed orbits with 
recurrence times T < Tmax is known in the energy interval -Emin < E < -Emax- The 
semiclassical eigenenergies and transition matrix elements are now obtained, in the 
second step, by adjusting the semiclassical signal to its quantum analogue (jTTj) 
with the {En, dn} being free adjustable parameters. The technical details to solve this 
nonlinear fit problem are given in Ref. [45]. 

The required signal length T^ax to achieve convergence of the harmonic inversion 
procedure depends on the mean level spacing g{E) in the energy range [£"111111, -^max] and 
reads Tmax > 4:7ig{E) [25]. The efficiency of the quantization method can be improved 
by using a cross-correlated semiclassical recurrence signal [25,26,46]. The idea is to use 
a set of L independent initial states I^Pj) and to construct the cross- correlated response 
function 

gf;^iE) = --{^,\DG+D\^,) ; z, j = 1, 2, . . . , L . (39) 



Application of the windowed Fourier transform as in Eqs. (jHTj) and ()38|) yields the 
quantum L x L cross-correlated time signal 

CT = -T.b^nb.ne-'''"', (40) 

J TT ^ — ' 

n 

with bin = (^il-Dl^n), and its semiclassical analogue 

Cr^{t) = — Ao,,(£)e^^-(^)e-^^*d£ . (41) 

•J Ejrnin CO 

For the various initial states the amplitudes Aco,ij{E) in equation (PT|) differ by 
the use of various angular functions yi^'d) in the closed orbit amplitudes ((Zj) [26]. The 
semiclassical cross-correlated time signal (PT|) can be adjusted to its quantum form (pUj) 
with the {En, bm} being the adjustable parameters by an extension of the harmonic 
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inversion method to cross-correlated time signals [25,47,48]. The idea is to identify the 
cross-correlated recurrence function C|j^(t), which is known on an equidistant time grid 
t = riT, with the cross-correlated time signal Cij{nT) = {^i\exp{—inTHea)\^j) of an 
effective Hamiltonian f^efr with the (not explicitly known) states |$j) and The 
operator Hes with eigenvalues En is complex symmetric, and {x\y) denotes a complex 
symmetric (not Hermitian) inner product. In an appropriate basis set the problem of 
extracting the {En, hn} can be reformulated as a generalized eigenvalue problem where 
all matrix elements can be expressed in terms of the time signal Cjj(nr). 

The advantage of using the cross-correlation approach can be understood based 
on the argument that the total amount of independent information contained in the 
L X L signal is L(L + 1) multiplied by the length of the signal, while the total number 
of unknowns {En,bin} is {L + 1) times the total number of poles En- Therefore the 
informational content of the Lx L signal per unknown parameter is increased, compared 
to the case of equation (jHEl), by a factor of L. This means that the required signal length 
Tmax = 4:'7ig{E) for a one-dimensional recurrence signal is reduced by about a factor of 
L for an L X L cross-correlated recurrence signal. 

In Sec. 14.21 we investigate dipole transitions from the initial state = \2p0) with 
light polarized parallel to the magnetic field axis to final states with magnetic quantum 
number m = 0. For this transition the angular function in equation (|7j) reads [25] 

yi{^) = (4 cos^ . (42) 

V 27r 

Results are obtained by harmonic inversion of a one-dimensional and a 2 x 2 cross- 
correlated recurrence signal. For the construction of the 2x2 cross-correlated signal we 
use for simplicity as a second transition formally an outgoing s-wave, viz. D\^2) oc yoo; 
and, thus, 3^2 ("i?) = 1- 

J^-.2. Results and discussion 

For the semiclassical photoabsorption spectra we calculated all closed orbits with 
recurrence times T < T^ax = 73.5 in the scaled energy range E e [—1, —0.5] on an 
equidistant energy grid with step width AE = 10~^. The closed orbits have been used 
to construct low-resolution photoabsorption spectra by superimposing the semiclassical 
contributions of the isolated orbits or, close to bifurcations, the uniform approximations. 
High-resolution photoabsorption spectra with individual semiclassical eigenenergies and 
transition matrix elements are obtained by the harmonic inversion method as explained 
in Sec. Wl\ 

Figure IHl presents the spectra for the photo excitation of the initial state |2p0) of the 
hydrogen atom in a magnetic field at field strength B = 11.75 T (7 = 5 x 10~^ a.u.) with 
light polarized parallel to the magnetic field axis. The semiclassical spectra in figure IHl^a) 
and (b) have been obtained by harmonic inversion of a one- dimensional recurrence signal 
and a 2 X 2 cross-correlated signal, respectively. In general, both semiclassical spectra are 
in good agreement with the quantum mechanical result. A detailed comparison shows 
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Figure 8. Semiclassical (SC) and quantum (QM) photoabsorption spectra of the 
hydrogen atom in a magnetic field at field strength B = 11.75 T (7 = 5 x 10~^ a.u.). 

Transitions from the initial state \2p0) with light polarized parallel to the magnetic field 
axis, (a) Semiclassical spectrum obtained by harmonic inversion of a one-dimensional 
recurrence signal, (b) Semiclassical spectrum obtained by harmonic inversion of a 2 x 2 
cross-correlated recurrence signal. 
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Table 1. Selected quantum and scmiclassical eigenenergies and transition matrix 
elements of the spectra shown in figure |H1 The indices s and c refer to semiclassical 
data obtained by harmonic inversion of a single (one-dimensional) and a 2 x 2 cross- 
correlated recurrence signal, respectively. 
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-1.24498 
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0.9813 
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that the cross-correlation technique improves the quahty of the semiclassical spectra. 
Some nearly degenerate states are not resolved with the one-dimensional signal but are 
well reproduced with the cross-correlation technique. Furthermore, the semiclassical 
and quantum transition matrix elements show better agreement in figure |Hfb) than in 
figure IHl^a). The eigenenergies and transition matrix elements of selected states are given 
in tabled The transition to the state at energy E = —1.22796 x 10^^ could not be 
resolved by harmonic inversion of the single recurrence signal (indicated by index s) 
but is well resolved when using the cross-correlated signal (see values with index c in 
tabled). 

In figure ini the magnetic field strength is reduced to i? = 4.7T (7 = 2 x 10^^ a. u.) 
which can be achieved easily in experiments. The same closed orbit data as in figure |H1 
have been used for the semiclassical calculations. Similar as in figure |H1 the semiclassical 
and quantum spectra in general are in good agreement, with the cross-correlation 
technique being even more reliable than the harmonic inversion of the one- dimensional 
recurrence signal. However, the convergence of the semiclassical spectra at B = 4.7 T is 
less perfect than at the higher magnetic field strength B = 11.75 T. The reason becomes 
evident from the scaling properties of the condition T^ax > 4:TTg{E) on the required signal 
length. If the scaled energy E is kept constant and the magnetic field strength 7 is varied 
the mean density of states scales {E = E72/3) = ~g^^^{E = E)7-^/3 whereas the 
recurrence time scales as T = T/7. Thus, in scaled units the required signal length reads 
Tmax > 47r^7~^/^, and becomes larger as the field strength is decreased. To improve the 
convergence properties of the harmonic inversion procedure closed orbits with longer 
periods are required, i.e, the cut-off limit Tmax for the scaled recurrence time must be 
increased. 
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Figure 9. Same as figure |H1 but for spectra at laboratory magnetic field strength 
B = 4.7T (7 ^ 2 X 10-5 a.u.). 
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5. Conclusion 

Almost a century after the postulation of Bohr's quantization rules for the hydrogen 
atom and a decade and a half after the emergence of closed orbit theory [10, 11] 
we have succeeded in calculating semiclassically from first principles high-resolution 
photoabsorption spectra of the diamagnetic hydrogen atom in the transition regime 
to chaos. The necessary tools, viz. closed orbit theory, uniform approximations at 
bifurcations, and the harmonic inversion method, although being known separately, 
have been combined for the first time to obtain individual semiclassical eigenenergies 
and transition matrix elements in that regime. 

The various steps can be summarized as follows: We have calculated all closed 
orbits in the energy range —1 < E < —0.5 with recurrence times T/27r < 12. The 
rotator orbits R'^ are created in a sequence of two bifurcations, viz. a pitchfork and a 
tangent bifurcation from the orbit perpendicular to the magnetic field axis. This rather 
complicated scenario is described by the normal form of the co dimension- 2 symmetric 
butterfiy catastrophe. The vibrator orbits are created in pitchfork bifurcations of the 
parallel orbit. Some of the symmetric vibrators and rotators undergo further pitchfork 
bifurcations, where pairs of asymmetric orbits are created. The pitchfork bifurcations 
are described by the normal form of the symmetric cusp catastrophe. For all bifurcations 
of closed orbits in the selected range of energies and recurrence times the uniform 
approximations have been constructed, which remove the divergences of the isolated 
orbit contributions. The contributions of the isolated closed orbits and the uniform 
approximations around bifurcations have been superimposed to obtain semiclassical 
low-resolution photoabsorption spectra, and, via a windowed Fourier transform, the 
semiclassical time signal C^^{t). The harmonic inversion method applied to that signal 
finally yields the high-resolution spectra with individual semiclassical eigenenergies and 
transition matrix elements. The method has been augmented by the cross-correlation 
technique to optimize its convergence properties and thus to further improve the quality 
of the results. Spectra have been obtained at magnetic field strengths B = 11.75 T and 
B = 4.7 T in the energy region —7^/^ < E < — 0.57^/^. The semiclassical spectra, 
especially those obtained with a 2 x 2 cross-correlated recurrence signal, show excellent 
agreement with the exact quantum spectra. 

The semiclassical calculations can, in principle, be extended to higher energies 
deep into the classically chaotic region of the diamagnetic hydrogen atom. However, 
this means that additional types of bifurcations and catastrophes must be considered 
for the construction of the uniform approximations, and, even worse, the numerical 
effort increases drastically due to the exponential proliferation of closed orbits in the 
chaotic regime. Clearly the objective of this paper was not to present a semiclassical 
method which is computationally more efficient than exact quantum computations (in 
fact, the opposite is true). Rather, the results are of fundamental importance as regards 
the development, understanding, and practical applications of semiclassical theories. 
These theories have already been successful in the hmiting cases of integrable and purely 
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hyperbolic chaotic systems. We have now closed the gap for systems with mixed regular- 
chaotic dynamics. 
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